Coherent Perfect Absorbers: Time-reversed Lasers 
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We show that an arbitrary body or aggregate can be made perfectly absorbing at discrete frequen- 
cies if a precise amount of dissipation is added under specific conditions of coherent monochromatic 
illumination. This effect arises from the interaction of optical absorption and wave interference, and 
corresponds to moving a zero of the elastic S-matrix onto the real wavevector axis. It is thus the 
time-reversed process of lasing at threshold. The effect is demonstrated in a simple Si slab geometry 
illuminated in the 500 — 900nm range. Coherent perfect absorbers are novel linear optical elements, 
absorptive interferometers, which may be useful for controlled optical energy transfer. 
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A laser is a physical system which, when subjected 
to an energy flux (pump), self-organizes at a threshold 
value of the pump to produce narrow-band coherent elec- 
tromagnetic radiation. In the absence of inhomogeneous 
broadening and quantum fluctuations, this radiation has 
zero linewidth. Above the first lasing threshold, lasers are 
non-linear systems, but at the first threshold they satisfy 
a linear wave equation with a negative (amplifying) imag- 
inary part of the refractive index, generated by the pop- 
ulation inversion due to the pump pQ. In conventional 
lasers, the gain medium is confined in relatively high 
quality-factor (high-Q) resonators, and the lasing modes 
are closely related to passive-cavity modes. However, re- 
cent demonstrations of random lasers [2j [3] have shown 
that the lasing threshold can be reached and coherent las- 
ing obtained in resonators with no high-Q passive-cavity 
modes. It can be rigorously shown within semiclassical 
laser theory that the first lasing mode in any cavity is an 
eigenvector of the electromagnetic scattering matrix (S- 
matrix) with infinite eigenvalue; i.e. lasing occurs when 
a pole of the S-matrix is pulled "up" to the real axis by 
including gain as a negative imaginary part of the refrac- 
tive index @]. This viewpoint suggests the possibility of 
the time-reversed process of lasing at threshold 16] . A 
specific degree of dissipation is added to the medium, cor- 
responding to a positive imaginary refractive index equal 
in absolute value to that at the lasing threshold. The sys- 
tem is illuminated coherently and monochromatically by 
the time-reverse of the output of a lasing mode, and the 
incident radiation is perfectly absorbed. We refer to such 
an optical system as a coherent perfect absorber (CPA). 

Coherent perfect absorption is a general and robust 
phenomenon related to the analytic properties of the S- 
matrix. For simplicity, we consider scattering in one or 
two dimensions, for which the electric field (in the TM 
polarization) is a scalar obeying the Hclmholtz equation: 

[V 2 +n 2 (r)fc 2 ] 4>{f) = 0. (1) 

Here k — lu/c (a scalar), uj is the frequency, and n — 
n' + in" is the complex refractive index, with n" < for 
gain and n" > for absorption. There is an "external 
region" , extending from some radius r s to infinity, where 



n = 1 (the following discussion can be extended straight- 
forwardly to the case of an external dielectric no). The 
field in the external region is a combination of incoming 
waves V'™ and outgoing waves ip^: 

0(f) = [«m C M + An C Ut (*=)]> r>r.. (2) 

in 

The scattering channel amplitudes are related by 
E m - S mm'Wa m ' = Pm- With n" = (lossless 
medium), S(k) is unitary for real k. Continuing k 
into the complex plane, S(k) possesses a countably in- 
finite set of poles and zeros, symmetrically placed at 
{&m =?mT «7m}, where 7 m > 0. 

When we add gain or dissipation, the zeros and poles 
of S(k) flow in the complex k plane, but they cannot 
simply appear or disappear, because arg [det(S)] always 
winds by 2n around each zero and —2tt around each pole. 
(We ignore pair creation and annihilation of zeros and 
poles, which normally does not occur.) Typically, adding 
dissipation moves the zeros down, and when there is suf- 
ficient dissipation a zero can cross the real axis at some 
specific value of the wavevector, k m . We parameterize 
the dissipation by the imaginary part of the refractive 
index, Im{n m } = n'^ > 0, and refer to the discrete 
pairs {{k m , n^J} as the CPA zeros. Radiation incident at 
each k m can be completely absorbed if it corresponds to 
the specific eigenvector of the S-matrix having eigenvalue 
zero. For homogeneous bodies such purely-incoming, 
real-fc radiation modes are the complex conjugate of the 
constant-flux states first proposed by Kapur and Peierls 
[7], and recently introduced into semiclassical laser the- 
ory by Tiireci et. al. Uj. We wish to emphasize that 
the CPA zeros are distinct from absorption resonances 
of the atomic or molecular medium, which do not re- 
quire specific illumination conditions. The CPA process 
arises from the interplay of interference and absorption: 
in the presence of specific amounts of dissipation, there 
exist interference patterns that trap the incident radia- 
tion for an infinite time. Even small rates of single-pass 
absorption can lead to perfect absorption, and hence me- 
dia that normally do not absorb radiation well at certain 
frequencies can be made to do so, albeit within narrow 



2 



lm(ka) 




FIG. 1: (color online) Motion of exact S-matrix zeros (blue 
crosses) and poles (blue triangles) in the complex-fc plane, as 
dissipation increases from zero, for a one-dimensional slab of 
length a and uniform index n. Also shown are the zeros (red 
circles) and poles (red squares) from Equation the single 
R-matrix pole approximation. Inset: schematic of system. 



frequency bands. We will demonstrate this in more detail 
below, with a simple realizable CPA in Si (silicon). Fi- 
nally, just as absorption can be enhanced by interference, 
we show that it can also be reduced by illuminating the 
CPA with eigenvectors of the S-matrix with constructive 
interference for escape. 

Before moving to specific examples, we discuss a gen- 
eral framework for finding CPA zeros based on the "R- 
matrix" theory of Wigner and Eisenbud [5] ■ In this ap- 
proach, the S-matrix is given by 

S(k) = ~e 2lkr ° [I - ikRik)}- 1 [I + ikR(k)} , (3) 

where tq > r s is an arbitrarily-chosen boundary far from 
the origin, and R{k) (the R-matrix) takes the form 

oo 

Rmm>(k) = C^C- (4) 

aa' — 1 

F a a' — {k a k ) S aa i ik 7ao'! (5) 



where tj>™ e R is the a-th Wigner-Eisenbud basis 
function evaluated at r = tq and decomposed into 
the m-th channel [10], while k a e 1 is the eigen- 
value of that basis function. The dissipation matrix, 
laa' — J d d r Im (n 2 ){r) <p a (r) </> a '(r), is real and positive- 
definite. 

When Q > 1, each S-matrix zero/pole is determined 
by approximating the R-matrix by a single term a; in 
this case the zeros and poles of the S-matrix occur at 



(1 + ilaa) k 2 Tif a k - k\ 



0. 



(6) 



where <p a = X)m(^™) 2 > 0- Without dissipation (j aa — 
0), this implies that all zeros have positive imaginary 
parts, as already noted. Furthermore, zeros cross the real 
axis exactly when j aa = ip a /k a , at frequency k = k a . 
In Fig. [l] we show the exact pole motion, for the case 



lm(n) 




Increasing k 



-»-o-»o+o-» at »-o-»-o 



.3.57 3.58 3.5 



2 

Re (n) 



FIG. 2: (color online) Complex refractive indices for the 
uniform dielectric slab as a coherent perfect absorber, for 
ka = 664.7. Parity-even solutions are shown as hollow blue 
circles; parity-odd solutions are omitted for clarity. The 
green curve shows the refractive index of Si at different fre- 
quencies. Inset: parity-even solutions (hollow blue circles) 
and parity-odd solutions (filled red circles) in the region 
3.55 < Re(n) < 3.62; the green cross shows the index of 
Si at ka = 664.7, a = 100 /zm (i.e. A = 945.3 nm). 



of the two-port CPA (see inset and discussion below), 
finding excellent agreement with the single-pole R-matrix 
prediction, even though the cavity only has Q ~ 30. 

The simplest possible CPA is a single port reflector, 
similar to the "critically-coupled fiber-resonator" sys- 
tems widely studied in integrated optics [TlTTCS] . How- 
ever, important properties of the CPA are only revealed 
when there are multiple input ports and hence non-trivial 
eigenvectors. We therefore study a two-port case to illus- 
trate the concept fully. Consider a finite one-dimensional 
slab of thickness a and uniform refractive index n, with 
two input channels for each k corresponding to incident 
radiation from left and right. It is straightforward to cal- 
culate the 2x2 S-matrix for arbitrary complex index n, 
and to show that it has a zero eigenvalue when 
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When ka 3> 1, we can find an infinite number of discrete 
solutions of this equation, n„ 
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(8) 
(9) 



In Fig. [2j we plot the solutions for ka — 664.7. Note 
that we have here restated the CPA problem so as to 
find the {n„(fc)} which produces a zero of the S-matrix 
for a fixed k. This can be achieved by letting both the 
real and imaginary parts of n v vary; earlier we varied the 
imaginary part n" at fixed n', leading to zeros at differ- 
ent fc-points. This is a useful reformulation because it 
suggests how a CPA can be achieved in practice. A real 
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FIG. 3: (color online) Semi-log plot of normalized output in- 
tensities (i.e. \s\ 2 , where s are S-matrix eigenvalues) vs. the 
wavelength A = 2n/k, for a 100 /im Si slab. Solid lines show 
log 10 (|s| 2 ) for a parity-even (blue) or parity-odd (red) eigen- 
mode. The dashed line shows 2(|r| 2 + |t| 2 ), the total output 
intensity when the two input beams are incoherent. Inset: 
upper and lower bounds for |s| 2 over a wide range of a. The 
gray area shows the actual bounds, and the solid black line 
show the approximate lower bound from Equation (10 1. 



material has a frequency-dependent n(k), and when the 
material is fabricated and illuminated in a 2-port configu- 
ration, one may realize a CPA by scanning k and looking 
for coincidences, i.e. n(k) w n v (k) for some integers v. 

In an indirect bandgap semiconductor such as Si, n" 
increases smoothly from very small values as k increases 
past the bandgap, while n' changes little [14]. For fre- 
quencies near the bandgap, dn" jdk dn"/dk and 
dn' jdk <C dn' v /dk; thus it is possible to vary k so as to 
pass very close to several CPA zeros. Like a bound state, 
a CPA zero has no intrinsic width; the scattered inten- 
sity near a CPA zero is I m i n ~ co\n(k) — n v (k)\ 2 , where 
Co is of order unity (neglecting the spectral width of the 
incident radiation). Fig. [3] shows the S-matrix eigenvalue 
intensities for a slab of undoped Si, with a = 100 /mi. 
For this system, \n(k) — n v (k)\ 2 < 10~ 5 for the optimal 
(u, k). As we see, deep CPA minima are achieved in the 
neighborhood of A = 945nm, with a maximum intensity 
contrast ~ 50 dB, and more than ten substantial minima 
are visible in the range 938nm < A < 954nm. Due to 
the mirror symmetry of the uniform slab, the two eigen- 
modes are parity eigenstates. The ± in Eq. ^ deter- 
mines whether the perfect-absorber condition is satisfied 
by the symmetric or antisymmetric mode. The location 
of these minima is a-dependent and hence tunable within 
a given material; we can derive a tight lower bound for 
the S-matrix eigenvalue intensities, 
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which goes to zero at the n values given in ^fh and locates 
the interesting operating regions. For Si, the optimal 
wavelength occurs at around 750 nm for a = 10 /im, 
and around 1000 nm for a = 150 um; the exact value 



depends on n(k), which in turn depends on the doping. If 
a > 50/Ltm, the spacing between CPA zeros, given in (|8| , 
becomes small, and the material index always passes close 
to one or more of them without fine-tuning a, making the 
behavior shown in Fig. [3] robust. Simulations of CPAs of 
this type for other indirect bandgap semiconductors such 
as GaP show similarly good results. 

Perfect absorption occurs only under the correct coher- 
ent two-channel illumination conditions, corresponding 
to the eigenvector of the S-matrix with eigenvalue zero. 
The effect arises from a combination of interference and 
dissipation: the reflected part of the first incident beam 
interferes destructively with the transmitted part of the 
second incident beam, and vice versa, and therefore the 
radiation is trapped in an interference pattern within the 
slab and lost entirely to dissipation. One sees from Fig. [3] 
that the other, orthogonal eigenmode has maximal |s| 2 , 
and in fact the material absorbs significantly less energy 
than if it is incoherently illuminated from both sides. In 
this mode, the reflected part of each beam interferes con- 
structively with the transmitted part of the other beam, 
causing the radiation to escape the slab more quickly 
than in the incoherent case. 

To illustrate the role of interference in the CPA, we 
write down the transfer matrix, a 2 x 2 matrix with unit 
determinant: 



T(n,k) = - 



1 



(11) 



where r(n, k) and t{n, k) are the reflection and trans- 
mission amplitudes for a single wave of unit amplitude 
incident from either direction. Perfect absorption occurs 
when Tn — 0, i.e. when r 2 = t 2 . (Furthermore, we can 
show that the non-zero eigenvalue, corresponding to the 
minimally-absorbed mode, is 2t.) In order for the two 
output beams to interfere destructively in the manner 
described above, not only must they have equal intensi- 
ties (|r| 2 = |i| 2 ), they must also have the correct rela- 
tive phase. This is impossible to satisfy when n is real, 
as expected from energy conservation. In the T-matrix 
analysis, this property is based on the impossibility of 
establishing the correct interference pattern; when n is 
real, T* = T^ 1 , and consequently (111]) implies that even 
if we achieve the condition |r| 2 = ftp, there is always a 
±7r/2 phase difference between r and t, so r 2 ^ t 2 . 

The above analysis applies to any 2-port system with 
mirror symmetry in the index variation. It can be ex- 
tended to arbitrary 2-port systems, for which the general 
transfer matrix has the form 



T = 



- TtTr 

-rt 



1 



(12) 



where rt and r r are the reflection amplitudes for a wave 
of unit amplitude incident from either the left and right 
respectively, and t is the transmission amplitude for ei- 
ther direction of incidence. Then the perfect absorption 
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condition can be written as rgr r = t 2 , and again this 
cannot be satisfied if the refractive indices are all real. 

One interesting property of the mirror-symmetric sys- 
tem is that we can change the output intensity simply 
by changing the relative phase of the input beams. For a 
phase- modulated input, [l,±e I?i ], the output intensities 
are equal in each channel, with the form 



,2 0.6 



I = I sin 2 



(13) 



For the uniform slab, Iq = |(n 2 — l)/(^ 2 + 1)| 2 , where n 
is the special refractive index satisfying ^ . The mirror- 
symmetric CPA therefore functions as a compact absorb- 
ing interferometer. Unlike an ordinary interferometer, it 
does not deflect the input beams into another output 
channel, but causes it to be absorbed entirely within the 
material. The CPA may therefore be good for applica- 
tions in which optical energy needs to be delivered to a 
device in a precise pulse shape or sequence, which could 
be controlled by modulating the relative beam phase in 
Eq. (13). Once the CPA has absorbed energy, it must 
flow out in some manner, e.g. as heat or as fluores- 
cence, the latter being less likely to be useful or desir- 
able. Hence direct bandgap materials, such as GaAs, 
which make good lasers, make poor CPAs since much of 
the absorbed light will be re-radiated. In contrast, in- 
direct bandgap systems such as Si, which are difficult to 
make into lasers, retain the energy as heat or as electron- 
hole pairs, leading to an enhanced photovoltaic response 
as CPAs. We emphasize however that if the incident ra- 
diation is broadband, then the response averages out to 
approximately the incoherent result, due to the oscilla- 
tion of the absorption in Fig. [3] 

In applications it will often be useful to have a large 
contrast between the perfectly absorbed and incoherent 
or reduced absorption illumination conditions. Accord- 



ing to ( 13 1, the intensity contrast between the maximally 
and minimally absorbing eigenmodes in a uniform slab 
is limited by n. For the Si slab of Fig. [3j the contrast is 
~ 0.75. (In particular, I < 1 for all n, because some field 
is always present in the absorbing material.) It is possi- 
ble to increase the contrast with a non-uniform system, 
such as a periodic array of slabs (i.e., a distributed Bragg 
reflector or one-dimensional photonic crystal). One such 
system, consisting of Si and non-absorbing Si02 slabs, 
is shown in Fig. |4j It exhibits a much larger absorption 
contrast than the uniform slab: the S-matrix eigenval- 
ues at the 956 nm absorption resonance have intensities 
\s 1 \ 2 = 5x 10~ 4 (symmetric) and |s 2 | 2 = 0.988 (antisym- 
metric). The absorption in the antisymmetric mode is 
exceptionally low because the field inside the slab is con- 
centrated in the Si02 regions, which are non-absorbing. 

Finally, we note the CPA is probably easiest to real- 
ize in a waveguide geometry, where it is relatively easy 
to control individual mode amplitudes. In a free-space 
geometry, it may be challenging to generate the input 



- v v v | 

Wavelength (nm) 


f 

\ 


920 


930 


940 


950 


960 




FIG. 4: (color online) Plot of S-matrix eigenvalue intensities 
s| 2 , as a function of input wavelength, for a mirror-symmetric 
"Bragg CPA" consisting of 20 Si slabs of width 188.5 nm, 
separated by SiC>2 slabs of width 266 nm. The refractive 
index of Si02 (n = 1.46) is real to a good approximation in 
this wavelength range, so absorption only occurs in the Si. 
A large (~ 99%) absorption contrast can be observed at 956 
nm, close to one of the system's photonic band edges. Inset: 
the intensity of the two eigenmodes, as a function of x, close 
to the edge of the system, at 956 nm. 



corresponding to the desired S-matrix eigenvector, due 
to the proliferation of channels. 
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